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Abstract 

We have proposed a procedure to extract the probability for valence particle being out of the 
binding potential from the measured nuclear asymptotic normalization coefficients. With this 
procedure, available data regarding the nuclear halo candidates are systematically analyzed and a 
number of halo nuclei are comfirmed. Based on these results we have got a much relaxed condition 
for nuclear halo occurrence. Furthermore, we have presented the scaling laws for the dimensionless 
quantity < r 2 > /R 2 of nuclear halo in terms of the analytical expressions of the expectation value 
for the operator r 2 in a finite square-well potential. 



1. INTRODUCTION 



Nuclear halo is a threshold phenomenon As the binding energy becomes small, the 
wave function of valence particle extends more and more outward if the centrifugal barrier is 
small (I < 1). Eventually, this leads to the wave function penetrating substantially beyond 
the range of nuclear force as the binding energy approches zero, i.e., occurrence of nuclear 
halo. 

The occurrence conditions for nuclear halo have already been discussed in detail in the 
literature DDDQ . Jensen and Riisagerp, recently, proposed the necessary conditions, 

A 2 

B p < 270— keV for s-states, (1) 

Zj 

Z < 0.44A 4/3 for p-states, (2) 

B p < 2A~ 2/S MeV for both s- and p-states, (3) 

where B p is the binding energy of valence particle. Eq.Q overrules Eq.Q and becomes 
decisive for all A for s-states as they pointed out. Most of halo nuclei experimentally 
observerd as well as theoretically predicted are s-wave halo. Hence, Eq.(jnj) is the most 
important conditions for halo occurrence. 

Among the nuclear halo candidates, however, only the well-known s-wave halos in 11 Li, 
u Be and the p-wave halo in n Be with the B P A 2 / 3 values of 1.48 MeV, 2.49 MeV and 0.89 
MeV, respectively, are below or close to the limit set by Eq.©. As will be seen , the others 
have binding energies much larger than this requirement. Therefore, it might be necessary to 
make a systematic inquiry concerning the occurrence conditions and scaling laws of nuclear 
halo. To this end, we have first to select, if possible, a proper experimental observable which 
can be employed to characterize the spatial extension of valence particle density. Direct 
measurements of the nucleonic density distribution should provide a visualized picture of 
nuclear skin and halo stuctures, therefore will be very interesting. Over the years a large 
variety of experimental methods have been developed, using leptonic probs (as electrons, 
muons, etc.) for investigating nuclear charge distributions, and hadronic probs (as protons, 
a— particles, pions, etc.) for exploring the distributions of nuclear matter p. All these 
methods were applied successfully for many years for the study of stable nuclei. Nowadays, 



study on the structures of exotic nuclei with many possible methods has become a new 
and exciting field of research. Recently, proton-nucleus elastic scattering at intermediate 
energies was applied for obtaining accurate and detailed information on the size and radial 
shape of halo nuclei Instead of density distribution, the probability of a valence particle 
being out of the binding potential may be a suitable quantity to assess the degree of the 
spatial extension. In the following, we will present a procedure to extract the probability 
from experimental data in a more reliable way, and try to get the occurrence conditions and 
scaling laws of nuclear halo, which might be in line with experimental observations. In this 
paper, we will limit ourselves to an investigation on two-body nuclear halos. 



2. PROBABILITY OF VALENCE PARTICLE BEING OUT OF THE BINDING PO- 
TENTIAL 

The probability for valence particle being out of the binding potential radius R can be 
evaluated by, 

, ^r 2 cg l3 {r)dr 
r{K, r 0} a ) — oo , (4 ) 

Jo r Kij{ r ) dr 

where <p n ij represents the normalized sigle-particle radial wave function in the (nlj) bound 
state, r and clq are the radius and diffuseness parameters of the potential,and the binding 
potential radius R is a equivalent square- well potential radius which can be derived from the 
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We have explicitly 



measured mean-square core radius by R = ~ (< r > c +4) fm 
indicated the dependence of the probability P on R, r and a®. At the asymptotic distance 
(finij behaves like, 

W_„, i+1/2 (2fcr) 

Pnlj - kj " , ( 5 ) 



where W_ v ^ z+i/2 ( r ) is the Whittaker function, k = ^j2[iB p /h 2 is the wave number, and 
/i and 7] are the reduced mass and the Coulomb parameter for system (A c + N). In the 
neutron case where i] = 0, Whittaker function is related to the modified Bessel function 
Ki+i/ 2 {kr) as W j + i/ 2 {2kr) = a/ krj txK ; + i/ 2 (At) [9] . In Eq. ©, by is the single-particle 
asymptotic normalization coefficient (ANC) defining the amplitude of the single-particle 
wave function in the asymptotic region. The single-particle ANC and the single-particle 
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spectoscopic factor are related to the nuclear ANC,C^ cNl j, bvjiol]. 

ci m = h 3 {s^y\ (6) 

for virtual process A ^ A c + N. Here A, A c and N stand for the nucleus, its core nucleus 
and the valence particle, respectively. Nuclear ANC is an experimentally measurable and 
independent of parameters of the potential Q, till • 

The probability calculated with Eq.(j3J) is only theoretical value and potential-parameter 
dependent. Eq. © sets a restraint on the single-particle wave function 0„^(r).The ampli- 
tude of 4>nij( r ) at the asymptotic distances, fry, is fixed as long as the values of C^ cNl j and 
are given. However, the potential parameters ro and ao, hence the single-particle wave 
functions are not determined uniquely. We search the potential parameters in such a way 
that the quantity, 

*p= 1^ [ S ij ^rdj (r<) - C AcNlj - I (7) 

becomes minimum. This is because (f)nij{f) has the asymptotic behaviour specified by Eq. (J3J) 
and Eq. © when Xp turns into minimum. We extract the probability from the experimental 
data of C A cNl j and as follows. First, single-particle wave function is calculated with 
the Woods-Saxon potential. The depth of the potential is adjusted to reproduce the binding 
energy B p . The values of diffuseness parameter ao are choosen in the range of 0.50 — 0.70 
fm. For each fixed a , the radius parameter r is varied in a small step till the minimum 
in xf, is reached. In the calculation of Xpi the summation runs from = 6.0 fm to 40 
fm in step of 0.1 fm . We find that the value of Xp does not change as long as Rn > 6 
fm. Then, for each pair of (r , ao) the probability for the valence particle being out of the 
binding potential radius R is calculated with Eq. (jlj). Alternatively, one may take the 
radius, R ws = r A l J^ [1 + 5/3(na /(r Al^ 3 )) 2 ] 1 ^ 2 |4| , as binding potential radius instead of 
R. We find that the value of R and the values of R ws at the minima in x\ are nearly the 
same. Hence, the probabilities calculated with Eq. (JHJ) below are nearly identical whether 
R or R ws is used. Besides, the results are the same within the experimental error when 
changes 10%. Because the relative uncertainty of the experimental is usually less 
than 10%, therefore, the results below are reliable within the accuracy of ^ ( ; p) . Fi srure 1 
shows the xl as a function of the probability P for the 2s excited state in 13 C. In this 



calculation, C% Nlj = 1.84 ± 0.16/m" 1 / 2 H and = 0.95 Q are used. It can be 



seen 



from the figure that the minimum in xt is very deep with a very small dispersion in the 



probability P. This is actually due to the fact that Xp reaches its minimum value when 
Eq. (jHJl is satisfied. Therefore, in terms of the experimentally measured nuclear asymptotic 
normalization coefficients C^ cNl j, we can calculate the average value of the probability P 

by, 

P = ^2w i P(R,r 0i ,a 0i ), (8) 

with 

Wi = — . (9) 

E (xjj ( r oi, aw)) 

The summation runs over the minimum points in \ 2 for different do- It is worth to note that 
the probability obtained in this way is nearly parameter independent. This is an interesting 
and meaningful result within the reach of present experimental knowledge. 

By means of the above procedure, the probability for the valence particle being out of 
the binding potential have been calculated for a number of nuclei. The values listed in Table 
1 are the weighted average probability P extracted from different experimental sources. 
The results are plotted in Figs. 2(a) and (b) function of 2^B p R 2 /h 2 and B P A 2 / 3 , 
respectively. Riisager et al. suggest a criterion for quantitative assessment of halos, 

i.e., the valence particle has large probability, say I 50%, of being out of the nuclear binding 
potential radius. According to this criterion, the nuclei under considerations are all s-wave 
halo in the ground state (solid points) or in the excited states (open circles). We find from 
Fig. 2(b) that halo may be able to occure for 

B p A 2/3 < 10 MeV (10) 

which is much relexed than the one given by Eq. (J3J). 

We have calculated the probability for the valence neutron in 2s-state using Woods-Saxon 
potential with normal parameters r$ = 1.27 fm, ao = 0.67 fm, U$ = (50 — 32 JV ^ 4 Z ) MeV 
for the nuclei with N = Z and (N — Z)/A = 1/3, respectively. The binding potential radius 
R used here is also calculated by R 2 = § (< r 2 > c +4) fm 2 and < r 2 > c = (r c^4^ 3 ) 2 ,the r 0c 
is obtained by fitting to the experimental rms radii of light nuclei. The calculated results are 
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also plotted in Fig. 2 as the solid and dashed lines. Basically, they are in agreement with the 
experimental data. In order to examine the effects of larger diffuseness, we have plotted in 
Fig. 2(a) the probabilities for the Woods-Saxon potential with the same parameters as those 
used above except for doubling the diffuseness, = 1.34 fm. As shown by the dash-dotted 
and dotted lines, these calculations overpredict the experimental data. It means that using 
a very large diffuseness in potential may not correspond to the realistic situation for nuclei 
with weakly bound neutrons. 

The probability of a valence particle being out of the square-well potential is, 



X + H Co 2 



P = X + 2 (l - ^ 
X 2 + 3 X + 3V 6 

For 1=2, 



I-72) ' = !■ (I 2 ) 



±r)-0.4 (x-0), (13) 
so 



(14) 



(15) 

In the above equations, U sq is the depth of the square-well potential. The dash-double- 
dotted line in Fig. 2(a) illustrats the square-well potential predictions for 2s-state. We see 
that it underpredicts the experimentally extracted data, implying the important roles of the 
potential diffuseness on the probability being out of the binding potential radius. 

3. SCALING LAWS OF TWO-BODY NUCLEAR HALO 
In terms of nuclear ANC, we can extract the root-mean-square {rms) radius of the 
probability distribution for valence particle in the orbit (nlj). It can be written as the 



(X + 1)(X 3 + 6X 2 + 15X + 15) 1 



where 



x = R 



2fiB p 



h 2 ' 

and 



2fiU sq 
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contributions from the interior and asymptotic regions [11, 1 



< r > 



1/2- 



R 



1/2 



Sij r^ nlj (r)dr + (Ci cNlj Y r*W\ (16) 

Jo JR 

The first term in the equation is somehow parameter dependent, while the second term is 
not. Moreover, in the case of weakly bound nuclei, the second term gives more than 90% 
contribution to the value of the rms radius. Thus the error introduced by the parameters 
is small in the cases under consideration. The rms radii of the valence particle have been 



calculated in this way for the nuclei 11 



15 



Be 0,0,0,13, 12b 111 



13,14,15Q 



11 



15 



19 



19 



C 



2l\ . They are listed in Table H] along with other parameters. Based on the assumption of 



a core plus a valence neutron structure, recently, Ozawa et al. |15( applied a Glauber-model 
analysis for a few body system (GMFB), and deduced spectroscopic factors £y for some 
selected nuclei from the measured interaction cross sections 07. With their parameters of the 
binding potential, we calculate the single-particle wave function and obtain single-particle 
ANC bij in asymptotic region. Then, the nuclear ANC can be obtained from the deduced 
S'-factor and the single-particle ANC with Eq. (jUJ). The rms radii for 11 Be, 15 ' 19 C are 
evaluated by means of the GMFB analysis, and the results are also listed in Table |1J In 
order to check the above results, the rms radii of the valence particle wave functions in 
these three nuclei are extracted by subtracting the core contribution from the mean-square 
matter radius 



A 2 



<r2> =AA h <r2> --T h <r2> " < 17 > 

where A, A c and Ah are the total, core and valence particle mass numbers of the system, 
respectively. These rms radii are compared with the other data in Table HJ Except for 
15 C, the rms radii of halo obtained with these three methods are in agreement within the 
experimental errors. 



Hamamoto and Zhang .228 have deduced the expressions for the expectation value of the 
operator r 2 in a finite square-well potential. The terms with £q in denominator in thier 
expressions are negligible in magnitude as compared to the other terms for the case of 
X 2 < 2 which we are interested in. After omitting them, we get the following scaling laws, 



< r z > 
R 2 



1 _ 



X 



I = 0, 



(18) 
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1 r/ 1 _ f (x+l) 2 +x + 3 + A_) + (\ + 1 >' 2 1 . < \ 2 + 2 V ! 2 ) 



i? 2 X 2 + 3x + 3 



SA 3 " 2 X ; 2£ 2 



3Co 2 

(19) 



Keeping the largest term of the above eqations in the limit x ~ ¥ 0, we will arrive at the 
scaling laws in Ref. Q]. It should be kept in mind that the above laws depend on the 
quantum number n through £q. Here n is the node number of the radial wave function of 
valence particle. If halo is defined in terms of the requirement that the experimental value 
of probability P is greater than 50%, or approximately x 2 < 1-8 (see Fig. 2 (a)) ,we get the 
following conditions from Eqs. (JTSJl and ()19|) for nuclear halo occurrence, 



< r 2 > 



R 2 



< r 2 > 



> 1.5 for 2s states, (20) 



> 1.9 for lp states. (21) 



R 2 

Since the probability P is less than 40% for I = 2, we come to the same conclusion as 
Riisager et al Q] that halo is unlikely to occur for the particle in the d states. 

In Fig. El the experimental data of < r 2 > /R 2 are compared with our scaling law as 
well as the predictions of the single-particle model for the valence particle in 2s state. In 
order to have a better statistics, we adopt the weighted average value of < r 2 > instead of 
the individual experimental results. For the data without error, we assigned them to 10% 
of uncertainty for evaluating the weighted average. We see from the figure that the scaling 
law Eq. (|18|l can account for the available experimental data of halo candidates,though it 
is derived in a finite square-well potential. Several authors jsl 0, 13, |3| have put forward 
their scaling laws. Being of their / and/or n independent, we do not present them in the 
figure. 

4. SUMARRY 

In summary, we have proposed a procedure to extract the probability for valence particle 
being out of the binding potential from the measured nuclear asymptotic normalization 
coefficients. With this procedure, available data regarding the nuclear halo candidates are 
systematically analyzed and a number of halo nuclei are comfirmed. Based on these results 
we have got a much relaxed condition for nuclear halo formation as compared to Ref.jfJ. The 
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effect of potential deffuseness on the probability being out of the nuclear binding potential 
radius is also discussed. In terms of the analytical expressions of the expectation value 
for the operator r 2 in a finite square-well potential, we have presented the scaling laws for 
the dimensionless quantity < r 2 > / R 2 of nuclear halo, which can account for the available 
experimental data of halo candidates. 

Acknowledgments 

This work was supported by the National Natural Science Foundation of China under 
Grants No. 10075077, 10105016, 10275092 and the Major State Basic Research Development 
Programme under Grant No. G200007400. 



[1] P.G. Hansen and A.S. Jensen, Ann. Rev. Nucl. Part. Sci. 45, 591 (1995). 
[2] K. Riisager, A.S. Jensen and P. Moller, Nucl. Phys. A548, 393 (1992). 
[3] D.V. Fedorov, A.S. Jensen, K. Riisager, Phys. Lett. B 312, 1 (1993). 

[4] D.V. Fedorov, A.S. Jensen, K. Riisager, Phys. Rev. C49, 201 (1994); Phys. Rev. C 50, 2372 
(1994). 

[5] A.S. Jensen, K. Riisager, Phys. Lett. B 480, 39 (2000). 

[6] P. Egelhof, G. D. Alkhazov, M. N. Andronenko, A. Bauchet, A. V. Dobrovolsky, S. Fritz, 

G. E. Gavrilov, H. Geissel, C. Gross, A. V. Khanzadeev, G. A. Korolev, G. Kraus, A. A. 
Lobodenko, G. Munzenberg, M. Mutterer, S. R. Neumaier, T. Schafer, C. Scheidenberger, D. 
M. Seliverstov, N. A. Timofeev, A. A. Vorobyov, and V. I. Yatsoura, Eur. Phys. J. A 15, 27 
(2002). 

[7] G. D. Alkhazov, M. N. Andronenko, A. V. Dobrovolsky, P. Egelhof, G. E. Gavrilov, H. Geissel, 

H. Irnich, A. V. Khanzadeev, G. A. Korolev, A. A. Lobodenko, G. Munzenberg, M. Mutterer, 
S. R. Neumaier, F. Nickel, W. Schwab, D. M. Seliverstov, T. Suzuki, J. P. Theobald, N. A. 
Timofeev, A. A. Vorobyov, and V. I. Yatsoura, Phys. Rev. Lett. 78, 2313 (1997) 

[8] K. Riisager, D.V. Fedorov, A.S. Jensen, Europhys. Lett. 49, 547 (2000). 
[9] M. Abramowitz and LA. Stegun, Handbook of Mathematical Function (Publications, 9th 
printing, New York, 1970). 



9 



[10] A.M. Mukhamedzhanov, C.A. Gagliardi, R.E. Tribble, Phys. Rev. C 63, 024612 (2001). 

[11] Z.H. Liu et al, Phys. Rev. C 64, 034312 (2001). 

[12] L. Jarczyk et al, Phys. Rev. C 28, 700 (1983). 

[13] K. Riisager, Rev. Mod. Phys. 66, 1105 (1994). 

[14] F. Carstoiu et al, Phys. Rev. C 63, 054310 (2001). 

[15] A. Ozawa et al, Nucl. Phys. A 691, 599 (2001). 

[16] F. Ajzenberg-Selove, Nucl. Phys. A 506, 1 (1990). 

[17] S. Fortier et al, Phys. Lett. B 461, 22 (1999). 

[18] T. Aumann et al, Phys. Rev. Lett. 84, 35 (2000). 

[19] Z.H. Liu, Chin. Phys. Lett. 19, 1071 (2002). 

[20] F. Ajzenberg-Selove, Nucl. Phys. A 523, 1 (1991). 

[21] P. Banerjee and R. Shyam, Phys. Rev. C 61, 047301 (2000). 

[22] I. Hamamoto, X.Z. Zhang, Phys. Rev. C 58, 3388 (1998). 

[23] PG. Hansen, B. Jonson, Europhy. Lett. 4, 409 (1987). 

[24] A.S. Jensen, E. Garrido, K. Riisage, and D.V. Fedorov, RIKEN Review 39, 3 (2001). 



10 



TABLE I: Deduced nuclear ANC, probability, rms radii and < r 2 > /R 2 for the nuclear halo 
candidates. 

a Deduced 



Nucleus 




Bp C^Nij 
(keV) (fm- 1 / 2 ) 


probability 
(%) 


< r 2 >l/2 

(fm) 


< r 2 > 1/2 
(fm) 


< r 2 > av /R 2 


Ref. 


u Be 


1/2+ 


504 


0.81±0.05 




6.68±0.43 a 






[15] 












6.65±0.31 b 






[15] 








0.76±0.03 




6.23±0.25 














0.78 




6.40 






[16] 








0.81 




6.68 






[12] 








0.78 




6.44 






M 










78.6±7.5 




6.46±0.16 


2.60±0.13 




12 B 


1- 


749 


0.94±0.08 


68.5±11.7 


5.64±0.90 


5.64±0.90 


1.86±0.60 


M 




2 - 


1696 


1.34±0.12 


54.2±4.0 


4.01±0.61 


4.01±0.61 


1.10±0.31 


UU 


13 C 


1/2+ 


1857 


1.84±0.16 


55.5±9.6 


5.04±0.75 


5.04±0.75 


1.55±0.47 


[UJ 




O - 


1274 


1.54±0.09 


64.0±7.5 


5.78±0.36 


5.78±0.36 


1.97±0.26 


[12i 




1- 


2083 


1.84±0.11 


56.1±6.7 


4.57±0.30 


4.57±0.30 


1.34±0.18 


[12] 




1/2+ 


1218 


1.05±0.22 




3.65±0.82 a 






[15] 












4.59±1.02 b 






[15] 








1.40 




5.40 






[20] 








1.49±0.15 




5.86±0.60 






[12 










60.4±6.6 




5.15±0.34 


1.55±0.21 




19 C 




240 


0.57±0.19 




7.87±1.49 a 






[15] 












7.63±2.46 b 






[15] 








0.55±0.07 




7.07±0.50 






[21] 










80.8±19.9 




7.17±0.47 


2.58±0.34 





with GMFB method. 6 Calculated with Eq. [T7| 
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FIG. 1: Dependence of the Xp on the probability for valence particle being out of the binding 
potential P for the 2s excited state in 13 C. Symbols are connected by a line for each ao value to 
guide the eye. 
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FIG. 2: Probability for valence particle being out of the binding potential as a function of 
2/iBpR 2 /h 2 (a) and B p A 2 / 3 (b). The solid points and open circles represent the s-wave halos 
in the ground state and in the excited states, respectively. The lines show the predictions of the 
single-particle models with Woods-Saxon potentials and the square-well potential. The arrows in 
the panel (b) illustrate the up-limits of B p A 2 / 3 value set by Eq. and Eq. (|lUj) . 
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FIG. 3: Experimental data of < r 2 > /R 2 vs 2/iB p R 2 /fr? for the valence particle in the 2s state. 
The dash-double-dotted line is the scaling law of Eq. (|18|) . The other lines show the results of 
the single-particle model calculations with Woods-Saxon potentials for the nuclei with N=Z and 
I=(N-Z)/A=l/3, respectively. 
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